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SECTION  1 


INTRODUCTION  AND  SUMMARY 

This  report  describes  a  theoretical  research  program  to  Identify  the 
mechanisms  of  Iso-lndex  behavior  In  materials,  l.e.,  the  crossover  wavelength 
of  the  ordinary  (no)  and  extraordinary  (ng)  refractive  Indices  In  certain 
uniaxial  crystals,  where  the  birefringence  (ng  ~  ng)  equals  zero. 

Research  of  this  nature  could  lead  to  an  understanding  and  Identification  of 
candidate  materials  suitable  for  use  In  very  wide  fleld-of-view,  very  narrow 
bandwidth,  optical  filters  for  strategic  laser  communications  In  the  blue/- 
green  spectral  region.  An  electro-optic  filter  designed  at  the  Hughes 
Research  Laboratories  (HRL)  relies  on  iso-lndex  uniaxial  crystals  which, 
because  of  dispersion  In  the  crystal  birefringence,  e:dilblt  a  zero  crossing  of 
birefringence  at  a  specific  wavelength  known  as  the  Iso-lndex  point.  The  Iso- 
lndex  phenomenon  has  been  observed  primarily  In  some  of  the  II-VI  compounds 
and  ternary  chalcopyrltes  of  the  type  which  fall  into 

hexagonal  (Including  trigonal)  or  tetragonal  crystal  classes. 

Our  work  in  this  area  proceeds  from  an  earlier  Navy  contract  (N00014-78- 
C-0201)  In  1978  to  Investigate  new  filter  approaches  to  meet  the  challenging 
requirements  of  satelllte-to-submarlne  communications.  That  study  led  to  our 
concept  of  the  Iso-lndex  coupled  wave  filter^  In  which  polarized  light  Is 
selectively  coupled  from  the  fast  to  the  slow  axes  In  those  crystals  exhibit¬ 
ing  a  zero  crossing  of  birefringence.  The  calculated  performance  of  this 
class  of  filters  showed  great  promise  for  meeting  blue/green  requirements. 

For  example,  a  filter  of  AgGaS2t  which  has  Its  Iso-lndex  point  at  4970  A, 
showed  a  calculated  bandwidth  of  0.3  A  over  a  full  2^  sr  FOV  for  a  1-cm  Inter¬ 
action  length. 

On  a  subsequent  program  (N00014-79-C-0491) ,  we  fabricated  test  devices 
and  demonstrated  the  narrow-bandwidth,  wlde-FOV  properties  offered  by  this 
technique.  Specifically,  our  AgGaS2  filters  with  3-mm  Interaction  length 
showed  a  0.94-A  bandwidth. 
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The  rapid  progress  In  this  area  suggested  that  It  was  timely  to  attempt 
to  determine  those  materials  chat  would  allow  this  technique  to  be  used  at  the 
candidate  laser  wavelengths  of  Interest  to  blue/green  programs.  The  general 
goals  of  this  study  effort  were 

*  To  understand  the  fundament  properties  of  blrefrlngent  crystal¬ 
line  materials  that  lead  to  Iso-lndex  behavior. 

*  To  develop  the  ability  to  predict  the  Iso-lndex  wavelength  and 
filter  bandwidth  for  new  Iso-lndex  materials. 

An  Important  specific  goal  of  this  program  was  to  provide  the  foundation  for 
Improved  Iso-lndex  filters  which  operate  at  shorter  wavelengths  than  currently 
possible  (wavelengths  from  4600  A  to  4800  A  are  desired)  and  exhibit  narrower 
bandwldths  than  are  currently  available  (bandwldths  as  narrow  as  0.1  A  are 
sought).  Such  filters  would  find  Immediate  use  In  strategic  laser 
communications . 

The  proposed  methodology  of  the  program  was  to  determine  energy  level 
diagrams  of  known  applicable  materials  of  Interest,  determine  polarization 
dependence  and  oscillator  strengths  of  the  bandgap  transitions;  then  calculate 
Iso-lndex  and  filter  characteristics  with  no  adjustable  parameters  for  known 
materials,  and  finally  predict  new  Iso-lndex  materials  and  their  potential  for 
narrowband  filter  applications. 

In  preliminary  studes,  we  conceived  and  pursued  a  phenomonologlcal  model 
which,  for  known  Iso-lndex  crystals,  correctly  predicted  the  location  of  the 
iso-lndex  point  and  the  strength  of  the  dispersion  at  that  location.  This 
model  recognized  that  the  large  dispersion  In  Index  of  refraction  near  the 
band  edge  Is  responsible  for  the  Iso-lndex  behavior.  Subsequently  we 
attempted  to  extend  this  work  to  Include  the  quantum  mechanical  nature  of  the 
band-edge  transitions  so  that  the  Iso-index  behavior  could  be  predicted  from 
first  principles  ~  a  feature  that  was  missing  In  the  earlier  theory.  It  soon 
became  apparent,  however,  that  It  did  not  appear  obvious  that  even  the  most 
highly  refined  first-principles  theory  would  In  fact  prove  universally  tract¬ 
able.  Also,  such  an  ambitious  approach  was  well  beyond  the  scope  and  capabili¬ 
ties  as  governed  by  time  and  funding  constraints.  We  therefore  decided  to 


10 


1 


proceed  by  refining  our  phenomenological  theory.  Using  several  different 
multiple-oscillator  models,  we  showed  that  It  Is  possible  to  predict  the  Iso- 
Index  wavelength  of  CdS,  for  example,  to  within  a  fractional  error  of  less 
than  O.IOZ.  In  addition,  we  demonstrated  the  capability  to  predict  the  ordi¬ 
nary  and  extraordinary  refractive  Indices  of  ZnS  over  a  broad  range  of  visible 
and  ultraviolet  wavelengths  to  within  a  fractional  error  of  less  than  0.50%. 
ZnS,  although  not  an  Iso-lndex  material.  Is  of  practical  Interest  to  this 
program. 

With  the  aim  of  developing  an  iso-lndex  material  having  a  cross-over 
wavelength  farther  Into  the  blue  than  that  of  CdS  (524.5  nm),  we  proposed 
to  study  the  system  Z]^_^Cd^S,  since  ZnS  has  a  greater  band  gap  than 
CdS  but  Is  nevertheless  Isomorphous  with  CdS.  The  reasoning  was  that  adding  a 
moderate  amount  of  zinc  would  shift  the  dispersion  curves  closer  to  the  blue 
while  still  preserving  the  Iso-lndex  characteristic.  Moreover,  there  has  been 
some  experimental  evidence  that  this  does  Indeed  occur  over  a  limited  composi¬ 
tional  range.  Our  theory  was  based  on  a  linear  Interpolation  between  the 
respective  dielectric  and  band  gap  parameters  of  ZnS  and  CdS.  At  the  present 
time,  the  theory  appears  to  require  further  refinement.  While  predicting  a 
monotonlc  decrease  In  wavelength  of  the  Iso-lndex  point  with  Increasing  zinc 
concentration  as  anticipated,  the  theory  does  not  agree  quantitatively  with 
the  limited  experimental  data  currently  available. 

In  Section  2  we  discuss  some  background  considerations  relevant  to  the 
theoretical  approaches.  Section  3  contains  the  development  of  our  Initial 
theoretical  work  based  on  a  single  classical  oscillator  model.  This  theory  Is 
further  refined  In  Section  4  In  which  a  more  physically  realistic  multiple 
oscillator  model  Is  developed,  leading  to  highly  accurate  values  for  the 
refractive  Index  of  CdS  and  ZnS  over  the  whole  range  of  visible  wavelengths. 

In  Section  5,  an  optional  model  based  on  the  highly  versatile  Fade  approxima¬ 
tion  to  series  solutions  for  the  refractive  Indexes  Is  put  forth,  yielding 
further  Insights  Into  both  the  strengths  and  limitations  of  the  theoretical 
models  as  applied  to  the  Zn-Cs-S  alloys. 
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SECTION  2 

BACKGROUND  THEORETICAL  CONSIDERATIONS 

In  1966,  C.H.  Henry  demonstrated  an  optical  filter  based  on  the  sign 
change  of  the  birefringence  that  occurs  near  the  band  edge  In  certain  uniaxial 
semiconductors.  He  showed  that  light  energy  can  be  coupled  between  the  ordi¬ 
nary  and  extraordinary  polarization  at  the  Isotropic  point  by  appropriately 
applied  strain  or  magnetic  field.  Thus,  when  placed  between  crossed  polar- 

Q 

Izers,  the  material  acts  as  a  narrowband  filter.  More  recently,  Laurentl  and 
coworkers  demonstrated  a  filter  tuning  capability  In  mixed  crystals  of  CdS  and 
CdSe.  At  HEIL,  we  have  demonstrated  how  wave  coupling  can  be  achieved  by  an 
applied  dc  electric  field,  and  we  discovered  a  highly  selective  filter  mate¬ 
rial,  AgGaS2,  for  the  blue/green  region.^ 

The  Iso-lndex  phenomenon  has  been  observed  primarily  In  II-VI  compounds 
and  ternary  chalcopyrltes**  of  the  type  addition,  a 

C  g 

few  examples  have  been  described.  Including  a  halide,  an  oxysalt,  and  an 
organic  crystal,  >  all  of  which  fall  into  the  tetragonal  or  hexagonal 
(Including  trigonal)  crystal  classes.  Known  materials  and  their  Iso-lndex 
wavelength  and  crystal  class  are  summarized  In  Table  1. 

An  examination  of  the  theoretical  explanation  for  the  Iso-lndex  phenom¬ 
enon  Introduces  several  concepts  Involving  relationships  to  band  gap,*^  com- 
presslonal  distortion,**  and  magnitude  of  the  birefringence.  Some  of  these 
concepts  are  discussed  below.  However,  while  a  unified  theoretical  picture  Is 
beginning  to  emerge,  none  exists  today  that  allows  one  to  predict  the  Iso- 
lndex  properties  of  new  crystals  or  existing  blrefrlngent  crystals  whose  Iso- 
lndex  properties  have  not  been  characterized.  One  observation  Is  the  decrease 
In  Iso-lndex  wavelength  with  Increased  band  gap.  We  plan  to  use  this  rela¬ 
tionship  to  tailor  the  blue/green  filter  currently  under  develoment  for  the 
4600  to  4800-A  region  (DARFA/ONR  Contract  No.  N00014-81-C-0783) .  Nonetheless, 
at  the  current  state  of  the  art  of  related  theory,  we  are  unable  to  explain 
and  predict:  (1)  whether  a  minimum  band  gap  exists  below  which  the  crossover 
does  not  take  place;  (2)  whether  there  Is,  therefore,  a  maximum  Iso-lndex 
wavelength  attainable;  and  (3)  what  the  required  basic  properties  are  to 
predict  new  materials. 
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The  existence  of  different  energy  gaps  In  blrefrlngent  crystals  for  light 
polarized  parallel  (II)  and  perpendicular  (1)  to  the  crystal's  c  axis  Is  well 
known.  It  has  become  apparent  that  this  property  leads,  under  certain 
circumstances,  to  the  Iso-lndex  effect.  Such  an  example  Is  shown  In  the  elec¬ 
troreflectance  data  of  Figure  1  for  AgGaSe2.  The  band-gap  transition  A  shown 
In  the  Inset  of  the  figure  occurs  at  the  lowest  energy  and  Is  predominantly 
excited  by  light  whose  electric  field  Is  parallel  to  the  c  axis.  Transition 
B,  which  Is  different  from  A  due  to  crystal  field  splitting.  Is  higher  In 
energy  and  predominantly  excited  by  light  polarized  perpendicular  to  the  c 
axis.  We  will  refer  to  these  two  energy  gaps  as  Eg  and  E^.  The  parallel 
direction  Is  associated  with  n^  (the  index  of  refraction  of  the  extraordi¬ 
nary  ray),  while  the  perpendicular  direction  Is  associated  with  (the 
Index  of  the  ordinary  ray). 

The  model  used  to  describe  the  Iso-lndex  behavior  Is  shown  In  Figure  2. 
The  figure  Includes  the  possibility  that  Eg  can  be  either  >Ex  or  <Ex  for 
both  positive  and  negative  uniaxial  crystals.  We  have  schematically  shown  both 
refractive  Indices  as  a  function  of  frequency  (v)  around  the  band  gap.  For 
simplicity,  higher  energy  transitions  have  been  Ignored.  An  increase  In  index 
occurs  because  of  resonance  enhancement  near  the  transitions.  The  birefrin¬ 
gence  An  =»  ng  -  Ug  is  also  shown  in  the  figure.  For  positive  uniaxial 
crystal  ( An  >  0)  at  frequencies  well  removed  from  the  band  gap.  It  Is  clear 
that  the  iso-index  point  (denoted  by  will  exist  only  for  dAn/dv  <  0. 

Thus,  a  necessary  condition  for  the  existence  of  an  Iso-lndex  material  Is  that 
the  sign  of  the  crystal  birefringence  must  be  opposite  to  the  sign  of  the 
slope  of  the  birefringence. 

One  group  of  materials  whose  birefringence  Is  well  characterized  near  the 
band  gap  Is  the  ternay  chalcopyrites.  These  materials  have  been  thoroughly 
studied  for  nonlinear  optical  devices.  Figure  3  shows  the  wavelength-depend¬ 
ent  birefringence  for  11  ternary  compounds.  The  figure  also  shows  the  data 
for  the  Il-VI  semiconductor  CdS  and  displays  three  of  the  four  possible  condi¬ 
tions  for  Figure  2.  The  following  observations  are  noted: 

•  All  the  II-IV-V2  crystals  (underlined  In  Figure  3)  are 
positive  uniaxial  with  positive  slope  and  therefore  do 
not  have  an  Iso-lndex  point. 
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Figure  1.  Eleccroreflectance  spectrum  of  AgGaSe?  (see  Ref.  4) 
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Figure  3.  Birefringence  of  ternary  chalcopyrites  (taken  from  Ref.  4). 

The  II-IV-V2  compounds  are  underlined.  Also  shown  is  the 
binary  II-VI  semiconductor  CdS. 
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*  The  I-III-VI2  class  contains  both  positive  and  negative 
uniaxial  crystals.  As  expected,  only  those  with  a  bire¬ 
fringence  whose  sign  is  opposite  to  that  of  the  slope  have 
Iso-lndex  points.  In  this  case,  they  are  all  negative  uni¬ 
axial.  Electroreflectance  data  verifies  that  En  >  Ej. 

for  these  crystals. 

*  An  example  of  a  positive  uniaxial  crystal  with  an  iso- 
lndex  point  is  the  binary  CdS.  Electroreflectance  data 
shows  that  Ej^  >  Eg  for  this  crystal. 

*  An  example  of  a  negative  uniaxial  crystal  with  no  iso-lndex 
point  is  CuInS2.  However,  low-temperature  reflectance  experi¬ 
ments  show  that  the  band  gap  transitions  are  polarization 
Independent  and  so  the  model  does  not  apply. 

The  model  described  by  Figure  2  could  be  useful  in  selecting  new 
materials  according  to  the  stated  necessary  conditions  and  provides  the  needed 
theoretical  foundation  for  studying  iso-lndex  materials.  The  theory  should 
Include  the  variations  in  oscillator  strength  that  are  graphically  implied  in 
Figure  2.  This  will  depend  ultimately  on  the  band  structure  of  the  semicon¬ 
ductors.  The  band  structure  is  affected  by  the  crystal  field  splitting  of  the 
valence  p  bonds,  departure  from  cubic  symmetry  (compresslonal  distortion),  and 
higher  lying  transitions.  An  Important,  and  perhaps  the  most  challenging, 
goal  is  to  gain  an  understanding  of  the  change  in  iso-lndex  behavior  with 
compositional  change  that  is  needed  to  tailor  iso-lndex  material  to  specific 
laser  wavelengths.  In  Figure  3,  for  example,  the  substitution  of  selenium  for 
sulfur  in  AgGaS2  shifts  the  AgGaSe2  iso-index  point  to  a  longer  wavelength,  as 
expected.  However,  the  entire  character  of  the  iso-index  material  CuGaS2  is 
changed  by  substituting  indium  for  gallium.  In  fact,  CuInS2  is  not  an  iso- 
lndex  material.  Both  are  negative  uniaxial  crystals,  but  in  this  case  the 
substitution  has  removed  the  polarization  dependence,  as  discussed  above. 
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SECTION  3 

SINGLE  OSCILLATOR  MODEL 

Preliminary  theoretical  calculations  of  the  performance  of  iso-index 
filter  materials  based  on  the  single  oscillator  model  described  in  Section  2 
yielded  surprisingly  good  results  in  predicting  the  iso-index  wavelength  and 
the  rate  of  change  of  birefringence  with  wavelength  in  the  vicinity  of  the 
iso-index  point. 

Ue  use  the  model  to  relate  the  slope  in  birefringence  to  optical  proper¬ 
ties  near  the  band  edge.  This  model  is  based  on  the  following  observations  of 
known  Iso-lndex  materials: 

*  There  is  a  slight  shift  in  the  band  edge  between  the  ordinary 
and  extraordinary  directions  of  the  crystal.  This  shift  is 
seen  in  the  reflection  data  of  Figure  1. 

*  The  index  of  refraction  away  from  the  band  edge  is  larger 

for  the  orientation  with  the  higher  energy  band  edge. 

*  The  iso-lndex  point  is  fairly  well  removed  from  the  band  edge. 

In  calculating  the  iso-lndex  point  and  its  related  slope,  the  last  observation 
allows  us  to  neglect  the  contribution  to  the  index  that  arises  from  absorp¬ 
tion.  The  first  two  lead  to  the  model  schematically  shown  in  Figure  4.  The 
refractive  index  of  the  higher  energy  orientation,  indicated  the  subscript 
1,  has  the  larger  index  of  refraction  away  from  the  band  edge.  The  parameters 
governing  the  iso-lndex  point  and  its  slope  are  -  ^2  (recall  the  dielectric 
constant  c  *  n  )  and  the  shift  in  the  band  edge. 

We  now  calculate  the  dependence  of  the  zero-crossing  and  its  related 
slope  on  the  governing  parameters.  The  dielectric  constant  of  the  materials  is 
taken  to  have  two  contributions.  Near  the  band  edge,  a  resonant  contribution 
arises  from  transitions  at  that  wavelength.  A  constant,  nonresonant  background 
term  is  used  to  account  for  all  remote  transitions  as  their  influence  is 
slowly  changing  over  the  range  of  Interest.  Thus, 
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where  1  refers  Co  Che  ordinary  or  exCraordlnary  dlrecClon,  is  Che  non- 
resonanc  concribuCion,  Is  che  suscepclblllcy  of  Che  reso- 
nanc  cranslclon  near  Che  band  edge,  and  >  c/\.  The  earlier  observaclon  chac 
Che  Iso-lndex  polnc  Is  fairly  well  removed  from  Che  band  edge  leads  Co  yec 
anocher  slmpllf IcaClon.  We  use  Che  Lorenczlan  llneshape  of  a  classical  oscll- 
lacor*  Co  describe  Che  resonanc  cranslclon.  This  llneshape  becomes 
Insenslclve  Co  Che  exacc  decalls  of  Che  llnewldch  and  falls  off  as  1/v  in  Che 
longer  wavelengch  region.  For  Che  cranslclon,  Che  suscepClblllCy 
Cakes  Che  form 


2Av 


H^i1 

-  V)  * 


(3.2) 


where  Che  Cerms  are  defined  In  Figure  4.  The  paramecer 
AVijAij  =  fjj  Is  Che  osclllacor  screngCh  of  Che 

cranslclon.  Summing  Equaclon  (3.1)  over  Che  resonanc  cransiclons  near 
gives 


(3.3) 


where  F,  ■  e  T  f 
i  n  ^ 


ij* 


We  now  calculaCe  Che  dlelecCrlc  consCanC  difference  beCween  Che  Cwo  crys- 
Cal  dlrecClons: 


Ae  > 


(3.4) 


For  slmpllclcy,  we  assume  chac  che  osclllacor  screngchs  In  Che  Cwo  dlrecClons 
are  comparable  and  Cake  Fj^  *  F2  ■  F.  The  success  of  chls  model  In  pre- 
dlcClng  che  zero-crossing  and  relaced  slopes  JusClfles  Chls  assumpclon,  ac 
lease  for  che  II-IV  semlconduccors  and  cemary  chalcopyrlCe  semlconducCors 
sCudled.  The  Iso-lndex  polnc  occurs  aC  Che  zero  of  EquaClon  (3.4),  which 
gives 


♦This  Is  an  approxlmaclon  applicable  when  exacC 

form  Is  shown  In  Equaclon  (4.1)  of  SecClon  4. 
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2F«\> 


e  -  e 
I  2 


2F6v 

s  - S- 

^1  ■  ^2 


! 

j 


(3.5) 


where  Wg  ■  (v^  +  '^)/2  la  the  average  band  gap,  and  6v  ■  -  V2  is 

the  different  In  band  gaps.  The  slope  Is  obtained  by  differentiating 
Equation  (3.4) 


d(Ae). 

dv 


o 


4f5'>  ('>  -  v  ) 

_ g  o' 


&  o 


(3.6) 

Then,  by  combining  Equations  (3.5)  and  (3.6)  and  again  neglecting  terms  of 
order  (d'^g)^,  we  get 


- 

dv 


''o> 


(3.7) 


Thus,  the  slope  Is  nearly  twice  the  difference  In  the  dielectric  constants 
divided  by  the  separation  of  the  zero-crossing  from  the  mean  band  gap. 

The  slope  of  the  birefringence  An  Is  related  to  that  of  Equation  (3.7) 
by 

Ae  »  +  n2)(nj^  -  n2)  *  2nAn 


Consequently,  Equations  (3.5)  and  (3.7)  become 


(^ 


o 


_g_4 


- 


(3.8) 


and 


df  An) 
dX 


2(nj  -  02) 

1  -  I 

o  g 


4cn 

F 


g  g  ■* 


1/2 


(3.9) 
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where  again  the  subscripts  o  and  g  refer  to  the  Iso-lndex  point  and  to  the 
band  gap,  respectively;  Is  the  average  band  gap;  and  5Xg  Is  the  shift 
In  band  gaps. 

We  now  proceed  to  calculate  the  zero-crossing  and  the  slope  of  birefring¬ 
ence  of  several  Iso-lndex  materials  and  to  compare  the  values  with  measured 
quantities.  Numerical  values  for  the  band-edge  shift  are  obtained  from 
reflection  data  such  as  those  shown  In  Figure  1.  The  crystal  birefringence  n^ 
-  n2  away  from  the  band  edge  Is  taken  from  previous  measurements  such  as  those 
In  Figure  3. 

The  measured  offset  of  the  Iso-lndex  point  from  the  band  edge  for  several 
II-VI  semiconductors  and  ternary  chalcopyrltes  Is  shown  In  Figure  5  along  with 
the  calculated  offset  from  Equation  (3.8).  The  agi 'between  theory  and 
experiment  In  Figure  5  Is  ~il5%,  which  we  feel  I*?  V  ;j.j;sd  considering  that 
there  Is  more  than  an  order  of  magnitude  varlaci^i.  >  *  ;d  gap  splitting 
and  birefringence  n^  -  n2  among  the  materials.  '^j^fncal  value  of  the 
constant  (F/cn)^/^  of  Equation  (3.8),  assumerf  o>;k.^ately  constant  for 
all  materials,  was  established  by  minimizing  t.t  :  deviation  of  the  calculated 
curve.  Thus,  the  factor  containing  the  oscillator  strength  appearing  In 
Equation  (3.3)  becomes 


cn  cn  -tZ- 
A  n 


(1.6  X  10"^)^  . 


(3.10) 


Evaluating  this  result  at  5000  A  as  an  example  gives  AX ‘A  •  160  A  for  n  »  2.5. 
Since  A  -  1,  this  result  shows  that  our  model,  which  neglects  the  llnewldth  AX 
of  the  transition.  Is  consistent  with  the  fact  that  (Xq  -  Xg)  »  AX. 

In  Figure  6,  the  measured  slope  Is  plotted  against  the  parametric  depend¬ 
ence  of  Equation  (3.9).  The  measured  slope  varies  by  a  factor  of  two  at  most 
from  the  calculated  values  over  the  range  of  three  orders  of  magnitude  of  the 
horizontal  axis.  We  believe  that  AgGaSe2  shows  the  most  disparity  because  It, 
unlike  the  other  three,  has  an  additional  band  edge  transition  not  accounted 
for  In  the  model  (transition  C  of  Figure  I). 

The  following  general  conclusions  can  be  drawn  from  this  work: 

*  Crystals  with  large  birefringence  are  the  best  candidates  for 
narrow-band  filters. 
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- - x/sx^  (W 

(n^  -  Oj) 


Xg(^) 

SXg(X» 

"l“"2 

X„-Xg(X) 

ZnO 

3637 

48 

0.013 

323 

AgGaS2 

4322 

422 

0.055 

648 

CdS 

4834 

38 

0.017 

411 

CuGaS2 

4973 

240 

0.006 

1341 

AgGaS«2 

6429 

686 

0.030 

1607 

Figure  5.  Measured  offset  of  the  Iso-lndex  point  from  band  edge  as 
function  of  material  parameters. 


Figure  6.  Measured  slope  at  Iso-lndex  point  as  function  of  material 


parameters. 


*  Crystals  such  as  Il-VI  semiconductors  with  smaller  band-gap 
shifts  are  preferable  If  the  crystal  birefringence  Is  suffi¬ 
ciently  large. 

*  For  comparable  birefringence  and  band-gap  shift,  the  shorter- 
wavelength  filters  offer  narrower  bandwldths. 

The  preceding  comparison  between  results  from  the  Iso-lndex  model  and 
measurements  gave  us  some  confidence  that  the  proposed  model  would  provide  the 
foundation  for  a  more  refined  theoretical  effort.  The  material  presented  In 
this  section  suggests  that  the  following  properties  of  candidate  uniaxial 
crystals  oust  be  known  In  order  to  predict  their  Iso-lndex  behavior: 

*  Energy  level  diagram  for  at  least  the  first  two  band-gap 
transitions 

*  Polarization  properties  of  these  transitions 

*  Oscillator  strengths  of  these  transitions 

*  Crystal  birefringence  (positive  or  negative  uniaxial) . 


SECTION  4 


MULTIPLE  OSCILLATOR  MODEL 

This  analysis  was  undertaken  to  Investigate  the  properties  of  the  indices 
of  CdS,  ZnS  and  their  alloys  Zn^.^Cs^S  to  predict  the  Iso-index  point 
and  the  slope  of  the  birefringence  as  a  function  of  the  composition  (>)■  A 
number  of  different  approaches  are  possible.  One  could  attempt  to  predict  the 
dielectric  constants  of  CdS  and  ZnS  from  a  band  theory  and  then  use  the  vir¬ 
tual  crystal  approximation  on  the  alloy  system  to  obtain  predictions  for  the 
alloys.  While  In  principle  this  approach  is  feasible,  a  review  of  the  suc¬ 
cesses  of  the  band  theory  In  predicting  the  optical  properties  of  these  mate¬ 
rials  Indicates  that  Che  success  of  this  approach  has  not  been  such  that  we 
could  expect  to  obtain  the  desired  degree  of  accuracy  for  the  dielectric 
constants. 

Q 

A  second  approach  Is  to  take  the  known  value  of  the  Indices  for  CdS  and 
ZnS  and  fit  them  to  a  physically  meaningful  form,  where  we  might  expect  the 
parameters  to  vary  linearly  between  the  two  cases  in  the  alloys.  This  is  the 
approach  that  we  followed. 

The  form  that  seems  most  reasonable  Is  the  sum  of  oscillators 

£(0.)  -  A(l)  +  y  - - r  .  (4.1) 

n»l  <«>(n)  -  <«>  +  A(o(n) 

In  this  equation,  the  A's  and  0J(n)*3  are  obtained  from  fitting  the  known  data.^ 
As  In  Section  3,  we  neglect  the  terms  d(u(n)  ln~  the  denominator  since  the  fre¬ 
quencies  of  interest  are  well  removed  from  the  band  edge.  The  symbology  of 
Equation  (4.1)  differs  from  that  of  Section  3  but  was  used  to  conform  to  that 
of  our  computer  prlnt-outs  displayed  In  the  accompanying  tables. 

This  approach  was  applied  using  one  and  two  oscillators.  (Trial  results 
using  more  than  two  oscillators  Indicated  no  improvement  In  matching  the 
formula  to  experimental  data.)  The  results  of  fitting  the  data  for  CdS  and 
ZnS  are  given  In  Tables  2  through  9.  As  can  be  seen,  the  fit  using  the  single 
oscillator  model  can  be  off  by  nearly  one  percent  while  the  calculation 
obtained  by  fitting  the  data  with  two  oscillators  Is  not  off  by  more  than  0.5% 
and  Is  typically  good  to  less  than  0.2%.  Hence,  the  two  oscillator  model 
gives  a  very  good  fit  to  the  data. 
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Table  2.  The  Result  of  Fitting  the  Data  for  the  Extraordinary  Ray  for 
CdS  to  the  Single  Oscillator  Model.  In  the  Table  EXP.  EP  and  TH.  EP  are 
Respectively  the  Experimental  and  Theoretical  Values  of  the 

Di' lectric  Constants. 


CdS  Extraordinary  Ray 
A(l)»  0.4454aD-«>l 
A(2)«  0.7ie83D-^01 
W(l)«  0.28«>d>8D-H:>l 


»ENGTH(A) 

FREQ. (aV) 

EXP.  EP 

TH.  EP 

7.  ERROR 

5120. 

2.421 

7.568 

7.504 

0.85 

5130. 

2.416 

7.524 

7.475 

0.66 

5140. 

2.411 

7.491 

7.446 

0.60 

5150. 

2.  407 

7.431 

7.418 

0.  17 

5160. 

2.402 

7.398 

7.391 

0.  10 

5170. 

2.397 

7. 366 

7.365 

0.  02 

5180. 

2.393 

7.322 

7.339 

-0.22 

5190. 

2.388 

7.301 

7.313 

-0.  17 

5200. 

2.384 

7.279 

7.289 

-0.  13 

5210. 

2.379 

7.258 

7.265 

-0.  10 

5220. 

2  .  o75 

7.231 

7.241 

-0.  14 

5230. 

2.370 

7.209 

7.218 

-0.  12 

5240. 

2.  365 

7. 182 

7.  195 

-0.  18 

5250. 

2.  361 

7.  156 

7.  173 

-0.25 

5275. 

2.350 

7. 102 

7.  120 

-0.25 

5300. 

2.339 

7.044 

7.070 

-0.37 

5op25  . 

2.328 

6.991 

7.022 

-0.45 

5350 • 

2.317 

6.954 

6.976 

-0.33 

5375. 

2.306 

6.906 

6.933 

-0.39 

5400. 

2.295 

6.875  ■ 

6.892 

-0.  25 

5425. 

2.285 

6.823 

6.852 

-0.44 

5450. 

2.274 

6.791 

6.815 

-0.  35 

5475. 

2.264 

6.760 

6.779 

-0.-2S 

5500. 

2.254 

6.724 

6.744 

-0.31 

5750. 

2.  156 

6.477 

6.467 

0.  15 

6000. 

2.066 

6.  305 

6.274 

0.49 

6250. 

1.983 

6.  170 

6.  132 

o;62 

6500. 

1 . 907 

6.  066 

6.024 

0.71 

6750. 

1 . 836 

5.983 

5.938 

0.75 

7000. 

1.771 

5.915 

5.869 

0.77 

7500. 

1.653 

5.803 

5.765 

0.66 

8000. 

1.549 

5.722 

5.690 

0.  55 

8500. 

1.458 

5.655 

5.635 

0 .  >^6 

9000. 

1.377 

5.607 

5.592 

0.28 

9500. 

1.305 

5 . 565 

5.558 

0.  12 

10000. 

1.239 

5.532 

'5.531 

0.  02 

10500. 

1.180 

5.504 

5.508 

-0.08 

11000. 

1.  127 

5.476 

5.489 

-0.25 

11500. 

1.078 

5.457 

5.473 

-0.30 

12000. 

1 . 033 

5.438 

5.460 

-0.40 

12500. 

0.992 

5.424 

5.448 

-0.44 

13000. 

0.953 

5.410 

5.438 

— 0 . 52 

13500. 

0.918 

5.396 

5.429 

-0.61 

14000. 

0.885 

5.387 

5.422 

-0.64 

30 


Table  3.  The  Result  of  Fitting  the  Data  for  the  Ordinary  Ray  for  CdS 
to  the  Single  Oscillator  Ilodel.  In  the  Table  EXP.  EF  and  TH.  EP  are 
Respectively  the  Experimental  and  Theoretical  Values  of  the 

Dielectric  Constants. 

CdS  Ordinary  Ray- 
Ad)-  0.4S217D-*-01 
A(2>«  0.57794D-*-Ol 
Wd)«  0.27S50D-^01 


WAVELENGTH (A) 

FREQ. <»V) 

EXP.  EP 

TH.  EP 

7.  ERROR 

5XS0. 

2.407 

7-524 

7.463  . 

0.78 

S160. 

2.402 

7.480 

7.432 

0.65 

3170. 

2.397 

7.437 

7.  399 

0.50 

3180. 

2.393 

7.388 

7.368 

0.26 

3190. 

2.388 

7.339 

7.337 

0.02 

3200. 

2.384 

7.301 

7.308 

-0.09 

3210. 

2.379 

7.290 

7.279 

0.  15 

3220. 

7.258 

7.251 

0.  10 

3220. 

7 . 220 

7.223 

-0 . 04 

3240. 

2.365 

7.  188 

7.  196 

-0.  12 

3250. 

2.361 

7.  150 

7.  170 

-0.28 

3275. 

2.350 

7.081 

7.  108 

-0.38 

5300. 

4.  •  W>W  7 

7.017 

7.049 

-0.  45 

5323. 

2.328 

6.959 

6.994 

—0 . 50 

5330. 

2,317 

6.906 

6.941 

-0 . 5 1 

5375. 

2. 306 

6.849 

6.892 

-0.63 

3400. 

2.295 

6.807 

6.845 

-0 . 56 

3423 . 

2.285 

6.770 

6.801 

-0.45 

5450. 

2.274 

6.729 

6.759 

-0.44 

3475. 

2.264 

6 . 69o 

6.718 

-0.39 

5300. 

2.254 

6. 656 

6.680 

—0 . 36 

3750. 

2.  156 

6.391 

6.380 

0.  16 

6000. 

2.066 

6.215 

6.  178 

0.59 

6230. 

1.983 

6.086 

6 .  Oww 

0.87 

6500 . 

1 . 907 

5.983 

5.  925 

0.  98 

6750. 

1 . 836 

5.890 

5.840 

0.86 

7000. 

1.771 

3.827 

5.  772 

0.  94 

7300. 

1.653 

5.712 

5.672 

0  .-70 

8000. 

1.549 

5.  636 

5.601 

0.62 

8500. 

1.458 

5.374 

5.543 

0.47 

9000. 

1.377 

5. 522 

5.  508 

0.  26 

9300. 

1.305 

5.480 

3.476 

0.07 

10000. 

1.239 

3.448 

5.451 

—0.06 

10500. 

1.  180 

5.420 

3 .  4 wO 

-0.  19 

1 1000. 

1.127 

5.401 

5.413 

-0.22 

1 1500. 

1.078 

3.382 

5.398 

-0.29 

12000. 

1 .  OO'O 

5.364 

5.386 

-0.  41 

12500. 

0.992 

5.345 

3.375 

-0 . 35 

13000. 

0.953 

5.331 

5.366 

-0.64 

13500. 

0.918 

5.318 

5.358 

-0.75 

14000. 

0.885 

5.308 

5 .  ^'5 1 

-0.79 
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Result  of  Fitting  the  Data  for  the  Extraordinary  Ray  for  CdS 
to  the  Two  Oscillator  Model.  In  the  Table  EXP.  EP  and  TH.  EP  are 
Respectively  the  Experimental  and  Theoretical  Values  of  the 

Dielectric  Constants. 

CdS  Extraordinary  Ray 
A(l)«  0.32103D+01 
A(2)«  0.77aiSD>00 
A  (3)-  0.2747H>-*-02 
W(l)»  0.2S675D+01 
W(2>«  0.37713D-^01 


WAVELENGTH (A) 

FREQ. (eV) 

EXP.  EP 

TH.  EP 

7.  ERROR 

3120. 

2.421 

7.568 

7.560 

0.11 

3130. 

2.416 

7.524 

7.318 

0.07 

3140. 

2.411 

7.491 

7.479 

0. 16 

3130. 

2.407 

7.431 

7.442 

-0. 15 

3160. 

2.402 

7.398 

7.407 

-0.  12 

3170. 

2.397 

7.366 

7.374 

-0.11 

3180. 

2.393 

7.322 

7.342 

-0.26 

3190. 

2.388 

7.301 

7.311 

-0. 14 

3200. 

2.384 

7.279 

7.282 

-0.03 

5210. 

2.379 

7.258 

7.253 

0.06 

3220. 

2.375 

7.231 

7.226 

0.06 

3230. 

2.370 

7.209 

7.200 

0. 12 

3240. 

2.363 

7. 182 

7.  175 

0. 10 

3230. 

2.361 

7. 156 

7.  151 

0.06 

’  3275. 

2.350 

7.102 

7.094 

0.11 

3300. 

2.339 

7.044 

7.042 

0.03 

3323. 

2.328 

6.991 

6.993 

-0.04 

3330. 

2.317 

6.954 

6.948 

0.08 

5375. 

2.306 

6.906 

6.906 

0.01 

3400. 

2.293 

6.875 

6.866 

0.  12 

5425. 

2.285 

6.823 

6.829 

-0..10 

3450. 

2.274 

6.791 

6.794 

-0.04 

5473. 

2.264 

6.760 

6.761 

-0.01 

3500. 

2.234 

6.724 

6.729 

-0.08 

3730. 

2.  156 

6.477 

6.479 

-0-03 

6000. 

2.066 

6.305 

6.305 

0.01 

6250. 

1.983 

6.170 

6.  173 

-0.04 

6500. 

1.907 

6.066 

6.069 

-0.04 

6730. 

1.836 

3.983 

5.984 

-0.01 

7000. 

1.771 

3.913 

3.913 

0.02 

7300. 

1 . 653 

3.803 

5.802 

0.01 

8000. 

1.349 

3.722 

3.720 

0 .  0^ 

8300. 

1.458 

5.655 

5.656 

-0.01 

9000. 

1.377 

3.607 

5.605 

0.05 

9500. 

1.303 

5.363 

5.564 

0.02 

10000. 

1.239 

5.532 

3.530 

0.04 

10300. 

1.180 

3.504 

3.301 

0.04 

11000. 

1.127 

3.476 

5.477 

—0 .  0^ 

• 

11300. 

1.078 

3.437 

5.457 

0.00 

12000. 

1.033 

5.438 

5.439 

-0.02 

12300. 

0.992 

3.424 

5.424 

0.00 

13000. 

0.953 

3.410 

5.411 

-0.01 

13500. 

0.918 

3.396 

5.399 

-0.05 

14000. 

0.885 

3.387 

5.388 

-0.03 

Table  5.  The  Result  of  Fitting  the  Data  for  the  Ordinary  Ray  for  CdS 
to  the  Two  Oscillator  Model.  In  the  Table  EXP.  EP  and  TH.  EP  are 
Respectively  the  Experimental  and  Theoretical  Values  of  the 
Dielectric  Constants. 

CdS  Ordinary  Ray 
A(l>»  0.25833D+01 
A<2)«  0. 10270D+01 
A(3)»  0.44S2SD>02 
W(l)-  0.2S567D+01 
W(2>»  0.42872D-^01 


WAVELENGTH ( A )  FREQ . 

(eV) 

EXP.  EP 

TH.  EP 

%  ERROR 

5150. 

2.  407 

7.524 

7.525 

-0.01 

5160. 

2.402 

7.480 

7.478 

0 . 03 

5170. 

2.397 

7.437 

7.434 

0.04 

5180. 

2.393 

7.388 

7.392 

-0.06 

5190. 

2.388 

7.339 

7.352 

-0.  19 

5200. 

2.384 

7.301 

7.315 

-0.  19 

5210. 

2.379 

7.290 

7.279 

0.  15 

5220. 

2.375 

7.258 

7.245 

0.  18 

5230. 

2.370 

7.220 

7.212 

0.  11 

5240. 

2.365 

7.  188 

7.  181 

0.  10 

5250. 

2.361 

7.  150 

7.  151 

0.00 

5275. 

2-.  350 

7.081 

7.081 

0.00 

5300. 

2.339 

7.017 

7.018 

0.  00 

5325. 

2.328 

6.959 

6.960 

-0.01 

5350. 

2.317 

6.906 

6.907 

0.00 

5375. 

2.306 

6.849 

6.858 

-0.  13 

5400. 

2.295 

6.807 

6.812 

-0.08 

5425. 

2.285 

6.770 

6.770 

0.01 

5450. 

2.274 

6.729 

6.730 

-0.02 

5475. 

2.264 

6. 693 

6.693 

0.00 

5500. 

2.254 

6.  656 

6.658 

-0.02 

5750. 

2.  156 

6.391 

6.391 

0.00 

6000. 

2.066 

6.215 

6.212 

0.05 

6250. 

1.983 

6.086 

6.080 

0.09 

6500. 

1.907 

5.983 

5.978 

0^09 

6750. 

1.836 

5.890 

5.895 

-0.07 

7000. 

1.771 

5.827 

5.826 

0.  03 

7500. 

1.653 

5.712 

5.718 

-0.  10 

8000. 

1.549 

5 .  6^6 

5.637 

-0.02 

8500. 

1.458 

5.574 

5.574 

0.01 

9000. 

1.377 

5.522 

5.524 

-0.03 

9500. 

1.305 

5.480 

5 . 483 

-0.06 

10000. 

1.239 

5.448 

5.450 

-0.04 

10500. 

1.180 

5.420 

5.422 

-0.04 

11000. 

1.127 

5.401 

5.398 

0.05 

11500. 

1.078 

5.382 

5.378 

0.09 

12000. 

1.033 

5.364 

5.360 

0.07 

12500. 

0.992 

5.345 

5.345 

0.01 

13000. 

0.953 

5.331 

5. 331 

0.00 

13500. 

0.918 

5.318 

5.320 

-0.04 

14000. 

0.885 

5.308 

5.309 

-0.02 

Table  6.  The  Result  of  Fitting  the  Data  for  the  Extraordinary  Ray  for  ZnS 
to  the  Single  Oscillator  Model.  In  the  Table  EXF.  EP  and  TH.  EP  are 
Respectively  the  Experlaental  and  TheoretlcaJ.  Values  of  the 
Dielectric  Constants. 


ZnS  Extreordinary  Ray 
A<1)«  0.34226Df01 
A<2)»  0. 376‘»7D-*-02 
W(l)»  0. 463790*01 


WAVELENGTH (A) 

FREQ. (aV) 

EXP.  Es 

TH.  EP 

7.  ERROR 

3600. 

3.443 

7.339 

7.327 

0.  16 

3730. 

3.305 

6-970 

6.984 

-0.21 

4000. 

3.099 

6-574 

6.588 

-0.22 

4100. 

3.023 

6-472 

6.470 

0.03 

4200. 

2.931 

6-376 

6.368 

0.  13 

4250. 

2.916 

6.320 

6.321 

-6. 02 

4300. 

2.883 

6-275 

6.278 

-0.05 

4400. 

2.817 

6-190 

6.200 

-0.15 

4500. 

2.734 

6.  136 

6.  130 

0.09 

4600. 

2.695 

6-  066 

6.068 

—0 .  Oo 

4700. 

2. 6^7 

6.017 

6.013 

0.08 

4730. 

2.609 

5.998 

5.987 

0.  13 

4800. 

2.582 

5.  968 

5.962 

0.  10 

4900. 

2.530 

5.919 

5.917 

0.04 

5000. 

2. ‘479 

5.881 

5.876 

0.08 

5250. 

2.361 

5.794 

5.788 

0.  10 

5500. 

2.254 

5.722 

5.717 

0.09 

3750. 

2. 156 

S.  655 

5.658 

-0.06 

6000. 

2.066 

5.607 

5.609 

-0.03 

6250. 

1.983 

5.  560 

5.367 

-0.  13 

6500. 

1.907 

5.522 

3.532 

-0.  17 

6750. 

1.836 

5.490 

5.501 

-0.20 

7000. 

1.771 

5.462 

5.474 

—0 .2^ 

8000. 

1.549 

5.420 

5.395 

0.45 

9000. 

1.377 

5.359 

5.345 

0.27 

10000. 

1.239 

5.304 

5.310 

-0.  11 

12000. 

1.033 

5.262 

5.267 

-0.08 

14000. 

0.885 

5.235 

5.241 

-0.  12 
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Table  7.  The  Result  of  Fitting  the  Data  for  the  Ordinary  Ray  for  ZnS 
to  the  Single  Oscillator  Model.  In  the  Table  EXF.  EF  and  TH.  EF  are 
Respectively  the  Experimental  and  Theoretical  Values  of  Che 
Dielectric  Constants. 

ZnS  Ordinary  Ray 
Ad)-  0.34165D-*-01 
A(2)»  0.37321D+02 
Wd)»  0.46306D+O1 


WAVELENGTH (A) 

FREQ. (mV) 

EXP.  EP 

TH.  EP 

7.  ERROR 

3600. 

3.443 

7.317 

7.309  • 

0.  1 1 

37S0. 

.  o05 

6.954 

6.965 

—0. 16 

4000. 

3.099 

6.554 

6.568 

—0 . 23 

4100. 

3.023 

6.447 

6.450 

-0.05 

4200. 

2.951 

6.360 

6.347 

0 . 20 

4250. 

2.916 

6.  305 

6.301 

0.06 

4300. 

2.883 

6.260 

6.258 

0.  03 

4400. 

2.817 

6.  180 

6. 180 

0.01 

■  4500. 

2.754 

6.116 

6.  110 

0.09 

4600. 

2.695 

6.047 

6.048 

-0.  02 

4700. 

2.637 

5.993 

5.993 

0 . 00 

4750. 

2.609 

5.978 

5.967 

0.  19 

4300. 

2.582 

5.944 

5.943 

0.02 

4900. 

2.530 

5.895 

5.897 

-0 . 04 

5000. 

Z.479 

5.861 

5.356 

0 . 09 

5250. 

2.361 

5.770 

5.768 

0.02 

5500. 

2.254 

5.693 

5.697 

-0.07 

5750. 

2.  156 

5.641 

5.639 

0.04 

6000. 

2.066 

5.584 

5.589 

-0.  10 

6250. 

1.983 

5.541 

5.548 

-0.  12 

6500. 

1.907 

5.504 

5.512 

-0.  16 

6750. 

1 . 836 

5.471 

5.482 

-0.20 

7000 . 

1.771 

5.438 

5.455 

-0.31 

8000. 

1.549 

5.401 

5.376 

0.'i6 

9000. 

1.377 

5 .  wo6 

5.326 

0.  19 

10000. 

1.239 

5.295 

5.291 

0.  06 

12000. 

1 . 0-j-^ 

5.244 

5.248 

-0 , 08 

14000. 

0.885 

5.221 

5.223 

-0.  03 
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Table  8.  The  Result  of  Fitting  the  Data  for  the  Extraordinary  Ray  for 
ZnS  to  the  Two  Oscillator  Model.  In  the  Table  EXF.  EP  and  TH.  EP  are 
Respectively  the  ExperlBsatal  and  Theoretical  Values  of  the 
Dielectric  Constants. 

ZnS  Extraordinary  Ray 
A<1)«  0.33450D+01 
A  (2)  —0.  19316D-02 
A(3)«  0.40192D>02 
W<1)«  0.34352D-^01 
W(2>«  0.46914D+01 


WAVELENGTH (A> 

FREQ. <»V> 

EXP.  EP 

TH.  EP 

X  ERROR 

3600. 

3.443 

7.339 

7.339 

0.00 

3750. 

3.305 

6.970 

6.969 

0.01 

4000. 

3.099 

6.574 

6.583 

-0.  14 

4100. 

3.023 

6.472 

6.467 

0.07 

4200. 

2.951 

6.376 

6.366 

0. 15 

4250. 

2.916 

6.320 

6.321 

-0.01 

4300. 

2.883 

6.275 

6.278 

-0.05 

4400. 

2.817 

6.  190 

6.200 

-0.  16 

4500. 

2.754 

6.  136 

6.  131 

0.07 

4600. 

2.695 

6.066 

6.070 

-0.05 

4700. 

2.637 

6.017 

6.014 

0.05 

4750. 

2.609 

5.998 

5.989 

0. 15 

4800. 

2.582 

5.968 

5.964 

0.06 

4900. 

2.530 

5.919 

5.919 

0.00 

50«30. 

2.479 

5.881  ’ 

5.878 

0.04 

5250. 

2.361 

5.794 

5.790 

0.06 

5500. 

2.254 

5.722 

5.719 

0.05 

5750. 

2.  156 

5.655 

5.660 

-0.08 

6000. 

2.066 

5.607 

5.610 

-0.05 

6250. 

1.983 

5.560 

5.568 

-0.  14 

6500. 

1.907 

5.522 

5.532 

-0.18 

6750. 

1.836 

5.490 

5.501 

-0.21 

70«X>. 

1.771 

5.462 

5.474 

-0.23 

8000. 

1.549 

5.420 

5.394 

0.46 

9000. 

1.377 

5.359 

5.343 

0 .  -^0 

10000. 

1.239 

5.304 

5.308 

— <?.08 

12000. 

1.033 

5.262 

5.264 

-0.03 

14000. 

0.885 

5.235 

5. 238 

-0.07 
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Table  9.  The  Result  of  Fitting  the  Data  for  the  Ordinary  Ray  for  ZnS 
to  the  Two  Oscillator  Model.  In  the  Table  EXP.  EP  and  TH.  EP  are 
Respectively  the  Experimental  and  Theoretical  Values  of  the 
Dielectric  Constants. 

ZnS  Ordinary  Ray 
A<1)»  0.336310+01 
A(2)»  0.784490-03 
A(3>«  0.390290+02 
W(l>»  0.344310+01 
W(2>«  0.466810+01 


WAVELENGTH (A) 

FREQ. (cV) 

EXP.  EP 

TH.  EP 

%  ERROR 

3600. 

3.443 

7.317 

7.317 

0. 00 

3750. 

3.305 

6.954 

6.935 

-0.02 

4000. 

3.099 

6.554 

6.565 

-0.  18 

4100. 

3.023 

6.447 

6.448 

-0.02 

4200. 

2.951 

6.360 

6.347 

0.22 

4250. 

2.916 

6.305 

6.301 

0. 07 

4300. 

2.883 

6.260 

6.258 

0.03 

4400. 

2.817 

6.  180 

6. 180 

0. 00 

4500. 

2.754 

6.116 

6.111 

0. 08 

4600. 

2.693 

6.047 

6.049 

-0.04 

4700. 

2.637 

5.993 

5.994 

-0,02 

4750. 

2^609 

5.978 

5.968 

0.  16 

4800 . 

2.582 

5.944 

5.944 

0.  00 

4900. 

2.530 

5.895 

5.899 

—0. 06 

5000. 

2.  479 

5.861 

5.858 

0.06 

5250 . 

2.361 

5.770 

5.770 

0. 00 

5300. 

2.254 

3.693 

5.698 

-0, 10 

3730. 

2.136 

5.641 

5.640 

0.02 

6000. 

2.066 

5.584 

5.590 

-0. 12 

6250. 

1.983 

3.341 

5.549 

-0. 13 

6500. 

1.907 

5.504 

5 . 5 1 

-0.-17 

6750. 

1.836 

3.471 

5.482 

-0.20 

7000. 

1.771 

5.438 

5.455 

-0. 31 

8000. 

1.349 

5.401 

5.376 

0.46 

9000. 

1.377 

5.  336 

5.325 

0.-21 

10000. 

1 . 239 

3.295 

5.290 

0.  09 

12000. 

1.033 

5.244 

5.246 

-0.04 

14000. 

0.885 

5.221 

5.221 

0.00 

The  results  of  calculating  the  Iso-lndex  point  and  slope  of  birefringence 
using  a  linear  interpolation  of  the  parameters  are  given  in  Table  10  for  the 
two  oscillator  model  and,  for  completeness,  for  the  single  oscillator  model  in 
Table  11.  The  results  for  the  two  oscillator  model  show  that  there  are  ranges 
of  X  where  no  iso-lndex  point  occurs.  For  other  values  of  x,  two  iso-lndex 
points  are  found.  For  those  x  values  where  two  iso-index  points  are  found, 
the  one  with  the  smallest  slope  and  the  same  sign  for  the  slope  as  the  other 
iso-lndex  points  is  used.  On  comparing  the  results  for  two  oscillators  with 
those  for  a  single  oscillator,  one  can  see  that  the  results  are  very  different 
even  in  the  regions  where  both  models  give  Iso-lndex  points.  This  fact  mili¬ 
tates  to  some  degree  against  the  validity  of  this  approach. 

The  results  for  the  Iso-lndex  point  and  slope  for  the  two  oscillator 
model  are  plotted  In  Figure  7.  The  data  show  a  monotonlcally  increasing  value 
of  the  Iso-lndex  wavelength  with  Increasing  x.  Some  experimental  data  are 
available  for  the  alloys.  These  data  also  give  a  monotonlcally  Increasing 
value  for  the  iso-lndex  wavelength  with  Increslng  x.  The  data  is  approxi¬ 
mately  linear. 

Xjso(*)  -  3877  +  983x(A)  for  0.4  <  x  <  0.7  .  (4.2) 

Our  results  are  approximately  linear  over  the  same  range  but  we  obtain 

XjSqCx)  -  3460  +  1583x(A)  for  0.40  <  x  <  0.7  .  (4.3) 

Hence,  we  differ  by  about  400  A  in  the  constant  and  almost  a  factor  of  two  In 
slope. 

This  discrepancy  remains  an  unsolved  problem  at  this  point.  There  Is 
some  doubt  as  to  the  accuracy  of  the  experimental  data  given  in  Reference  3, 
but  we  have  no  way  of  quantifying  this.  With  respect  to  our  theoretical 
approach,  we  feel  that  pushing  the  empiricism  as  far  as  we  have  may  not.  In 
fact,  be  warranted.  Reliable  experiments  are  needed  to  yield  valid  Informa¬ 
tion  on  the  variation  of  the  Iso-lndex  point  with  alloy  composition. 


Table  10.  The  Iso-Iadez  Point  and  Slope  aa  a  Function  of  Coapoeition  for  the 
Two  Oscillator  Model.  WISO  is  the  Iso-lndex  Frequency  Measured  in  Electron 
Volts;  LISO  is  the  Iso-Index  Wavelength  Measured  in  A.  ODn/Dw  and  OOn/DL  are 
the  Derivatives  of  the  Difference  in  the  Indices  at  the  Iso-Index  Point  with 
Respect  to  Frequency  and  Wavelength,  Respectively.  The  Parameters  used  are 
at  the  Top  of  the  Table.  The  Range  in  which  no  Roots  were  Found  is 

Measured  in  Electron  Volts. 


•  • 

CdSC 

CdSO 

ZnSE 

Zn80 

A(l>- 

A(2)* 

.4AS48e-»01 

.71883C«Oi 

.45217E'K>1 

.57744e-*01 

.34226E4-01 

.37647E+02 

.341658*01 

.573218*02 

U(t>- 

.28668C«01 

.27e50E*01 

.46574E-»01 

.46306E*01 

Zn<t-i«) 

Cd<K)S  On^ne^na 

- 

X 

MlSO<aV> 

LlSOlA) 

OOn/OH 

DDn/DL 

0.00 

3.416 

3164.4 

. 151S2E-01 

-. 187B7E-04 

.05 

3.544 

3442.8 

.40721B-02 

-.42172E-05 

.to 

3.184 

3887.1 

.SB484E-02 

-.474818-05 

.15 

2.ai7 

4400.4 

.37851E-02 

-.242168-05 

.20 

2.408 

5148.0 

.232038-02 

-. 10852E-O5 

.25 

1.410  • 

6440.4 

. 122B4e-02 

-.3613aE-06 

.30 

1.047 

11247.7 

.37562E-03 

-.36477E-07 

.  35 

NO  ROOTS  MCRE 

ROUND  IN  .400 

TO 

3.485 

.40 

NO  ROOTS  HERE 

ROUND  IN  .400 

TO 

3.842 

.45 

NO  ROOTS  MERE 

ROUND  IN  .400 

TO 

3.800 

.50 

NO  ROOTS  HERE 

ROUND  IN  .  400 

TO 

3.708 

.55 

NO  ROOTS  HERE 

ROUND  IN  .400 

TO 

3.615 

.60 

NO  ROOTS  HERE 

ROUND  IN  .400 

TO 

3.523 

.65 

NO  ROOTS  HERE 

ROUND  IN  .400 

TO 

3.431 

.70 

1.434 

S614.3 

. 182368-02 

-.30461E-O6 

.75 

l.SOO 

6885.2 

.543648-02 

-.142148-05 

.SO 

2.006 

6174.1 

.118628-01 

-.58506E-05 

.85 

2. 144 

5782.0 

.233288-01 

-.86488E-05 

.40 

2.242 

5524.6 

.442868-01 

-. 17452E-04 

.45 

'  2.313 

5354.8 

.845018-01 

-.36460E-04 

t.OO 

2.364 

5242.8 

. l*B44e*CO 

-.  754548-04 

39 


Table  11.  The  Iso-Iadex  Point  and  Slope  as  a  Function  of  Composition  for  the 
Single  Oscillator  Model.  WISO  Is  the  Iso-Index  Frequency  Measured  In 
Electron  Volts;  LISO  Is  the  Iso-Index  Wavelength  Measured  In  A.  DDn/Dw  and 
DOn/DL  are  the  Derivatives  of  the  Difference  In  the  Indices  at  the  Iso- 
Index  Point  with  Respect  to  Frequency  and  Wavelength,  Respectively.  The 
Parameters  used  are  at  the  Top  of  the  Table.  The  Range  In  which  no  Roots  were 

Found  Is  Measured  In  Electron  Volts. 


CdSf 

CdSO 

Zn86 

ZnSO 

A<1)* 

.32103e«01 

258336*01 

.334506*01 

.336316*01 

A(2>a 

.778186*00 

« 

102706*01 

.OOOt^l56*<>l 

.  784496-05 

A(3>- 

.274716*02 

• 

448256*02 

.401926*02 

. 390296*02 

»(!>■ 

.256756*01 

2SS67E*bl 

-34352E*<M 

.344316*01 

M(2>a 

.377136*01 

• 

428726*01 

.469146*01 

.466816*01 

Zn kl-M iCd (x >S  0nane-n« 

X 

MlSO(dV) 

LISO (A) 

ODn/DN 

DDn/DL 

0.00 

NO 

ROOTS  WERE 

FOUND 

IN  .VOO 

TO 

3.435 

.OS 

NO 

ROOTS  MERE 

FOUND 

IN  .V<Xt 

TO 

3.392 

.  10 

NO 

ROOTS  MERE 

FOUND 

IN  .VOO 

TO 

3.348 

.  IS 

NO 

ROOTS  MERE 

FOUND 

IN  .900 

TO 

3.305 

. 

.20 

NO 

ROOTS  MERE 

FOUND 

IN  .900 

TO 

3.262 

.25 

3.  181 

3896.3 

. 408736*00 

-.3337lE-*>3 

.2S 

3.203 

3870.1 

-.20198E*»>1 

. 167156-02 

.  30 

3.126 

3964.5 

.516826*00 

-.407586-03 

.35 

3.074 

4032.2 

.55340E*00 

-.42188E-03 

.40 

3.  022 

4101.4 

.567136*00 

-.417896-03 

.40 

3.087 

4015.5 

-.49425E*03 

. 379946*00 

.45 

2.V70 

4172.9 

.565726*00 

-.402696-03 

.50 

2.  V19 

4246.9 

.555976*00 

-.382086-03 

.55 

2.867 

4323. 9 

.53927E*00 

-. 357526-03 

.60 

2.814 

4404.3 

.517136*00 

-.330446-03 

.65 

2.762 

4488.4 

.49135E*00 

-.302316-03 

.70 

2.708 

4577.1 

.460026*00 

-.272176-03 

.75 

2.654 

4670.5 

.42698E*00 

-.242616-03 

.80 

2.5VV 

4769.8 

.391106*00 

-.213086-03 

.85 

2.542 

4875.8 

.353166*00 

-. 184136-03 

.VO 

2.484 

4989.6 

.314216*00 

-.  156446-03 

.V5 

2.424 

5113.1 

.274496*00 

-. 130146-03 

1.00 

2.361 

5248.9 

.234166*00 

-. 105346-03 

40 
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SECTION  5 

THE  FADE  APPROXIMATION 


The  Fade  approximation  Is  a  very  useful  and  powerful  method  for 
representing  any  analytic  function  (which  Is  more  traditionally  represented  by 
an  Infinite  series)  In  terms  of  a  quotient  of  two  series,  usually  of  much  more 
limited  order.  Since  the  real  part  of  the  refractive  Index  Is  an  analytic 
function  of  the  frequency  squared  (^^  ),  we  attempted  to  fit  the  four  Indices 
of  refraction  of  CdS  and  ZnS  tabu  ated  by  Blenewskl  and  Czyzak^  using  the  form 

N  . 

V  ,.2n 

Z  ^n“ 

nW  =  -  •  (5.1) 

1  +  ^  b/“ 

m=l 


The  method  used  to  fit  a  set  of  data  points 

^x^iy^Jj  1  •  1,2,...,  I 


by 


>  a  X, 
z  n  1 


n*0 


M 


1  + 


y  bx“ 

^  mi 


m»l 


(5.2) 


Is  derived  multiplying  each  side  by 


1  + 


yb  x" 
^  m  1 


;  k  -  0,  I,  2,  ...,  N  +  M 


and  summing  over  1.  This  yields  a  set  of  simultaneous  linear  equations  for 
the  Sq  and  hg,: 
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M 


2  2  *a  •  ^1^  -  1  '’mf  ° 


a*0 


m-l 


(5.3) 


which  may  be  solved  directly.  Thus  the  method  Is  far  more  efflclenct  than  a 
direct  nonlinear  least-square  error  determination  of  the  coefficients. 

We  also  attempted  to  fit  (Re  n)^  as  P((o)/Q(ia)  >  as  opposed  to  fitting 
Re  n  to  P(‘A^)/Q(<d^)i  but  the  results,  for  the  same  number  of  parameters,  were 
not  as  good. 

It  was  found  that  Increasing  N  and  M  In  Equation  (5.1)  beyond  2  gave  no 
appreciable  Improvement  In  the  fit.  The  rms  and  maximum  discrepancies  for  the 
N  *  2,  M  >  2  case, 

2  4 

a  +  a.m  +  a_m 

n(U)2_0 - 1^ - 2^ 

1  +  bj^m  +  b^w 


rms 

max 

ZnS  -  extr. 

0.0019 

0.005 

ZnS  -  ord. 

0.0019 

0.005 

CdS  -  extr. 

0,0012 

0.004 

CdS  -  ord. 

0.0010 

0.003 

The  lack  of  improvement  for  Increased  M  and  N  and  the  generally  "random" 
distribution  of  the  residuals  suggests  that  they  are  due  principally  to 
measurement  errors.  There  was,  however,  some  Indication  of  (narrow  band) 
Impurity  absorption  between  0.7  and  0.8  ta  In  the  ZnS. 

The  values  of  the  coefficients  that  fit  the  four  n(u^)  functions  were: 


a 

0 

a^Csac  ) 

4 

a^Csec  ) 

K  ,  2. 

bj^Csec  ) 

bjCaec^) 

ZoS-a 

2.2743 

-8.3963  X  lO"^^ 

7.3132  X  lO"^* 

-4.0365  X  10“^^ 

3.9840  X  10~®^ 

ZnS-o 

2.2716 

-9.3590  X  lO"^^ 

9.3463  X  10 

-4.4531  X  lO"^^ 

4.9478  X  lO"** 

CdS-e 

2.2921 

-1.7755  X  lO"^^ 

2.1257  X  lo"*^ 

-8.4068  X  lO"^^ 

1.3055  X  lO"®^ 

CdS-o 

2.2763 

-1.8735  X 

2.5311  X  lO"*^ 

-8.8744  X  lO"^^ 

1.5003  X  10 

with  ■  2‘^c./^. 
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The  ratio  of  a  fit  quadratic  In  <0^  Is  equivalent  to  a  "two  oscillator" 


model : 


where 


n((0) 


a  +  a,  +  a-W^ 

3  _o _ 1 _ i 

2  4 

1  + 


(5.5) 


The  values  of  e,  fi^2»  “l,2  thus  obtained  were: 


e 

f^(sec"^) 

Mj^(sec 

f2(sec 

ti)2(sec”^) 

ZnS-e 

1.8356 

2.1402  X  10^° 

6.5697 

X  10 

2.2631  X  10^^ 

7.6261 

inl5 
X  10 

ZnS-o 

1.8890 

-1.2677  X  10^^ 

6.5589 

X  10^^ 

3.823  X  10^^ 

6.8544 

X  10 

CdS-e 

1.6283 

5.4890  X  10^^ 

3.9680 

X  10^^ 

3.0601  X  10^^ 

6.9751 

X  10 

CdS-o 

1.6870 

1.5148  X  10^^ 

3.8920 

XI0I5 

2.4639  X  10^^ 

6.6335 

1015 
X  10 

The  negative  value  for  f^  of  ZnS-o  Is  puzzling,  since  It  has  no  obvious  physic 
meaning. 

Lacking  any  guidance  from  first  principles  on  how  to  Interpolate  the  data 
from  CdS  and  ZnS  to  estimate  the  refractive  properties  of  Znj^.^Cs^S,  we 
tried  four  approaches:  (1)  Interpolation  of  e,  fi^2»  “l,2^» 

(2)  Interpolation  of  a^^,  b^,  (3)  Interpolation  of  Re  n,  and  (4)  Interpola- 
tlon  of  (Re  n)  .  (By  "Interpolation"  we  mean  linear  Interpolation: 

Qx  ■  *QcdS  +  (1  "  *)  QznS*^  Interpolation  of  the  oscillator 

(method  (1)),  which  we  had  anticipated  would  be  a  plausible  approach,  gave 
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poor  results.  Specifically,  the  Iso-lndex  point  was  "sucked"  Into  the 
absorption  region  (near  t»^)  for  values  of  x  very  near  unity.  Interpolation  of 
the  coefficients  a^,  did  give  results  that  were  at  least  believable. 

The  Iso-lndex  point  decreased  In  wavelength  with  decreasing  x,  crossing  Into 
the  absorption  region,  however,  at  x  >  0.76,  X  =  4340  A.  An  Iso-lndex 
wavelength  of  4800  A,  for  example,  was  predicted  for  x  "  0.83.  (In  the 
approach  of  Section  4,  an  Iso-lndex  wavelength  of  4800  A  requires  x  "  0.80,  as 
Indicated  by  Figure  7.)  Interpolation  of  Re  n  and  (Re  n)^  gave  results 
similar  to  each  other,  with  the  slope  dX^^^/dx  again  positive,  but  much 
smaller  than  that  of  method  (2). 

The  results  obtained  \sy  method  (2)  are  given  In  Figure  8.  Included  In 
Figure  8(a)  for  clarification  are  the  two  lower  frequency  poles  of 
Equation  (5.4).  It  Is  seen  that  the  Iso-lndex  wavelength  crosses  Into  the 
absorption  region  near  x  *  0.76,  as  mentioned  above.  As  a  mathematical  conven¬ 
ience,  we  plotted  In  Figure  8(b),  the  slope  of  birefringence  multiplied  by  the 
wavelength  at  the  Iso-lndex  point.  The  sharp  Increase  In  this  slope  Is  to  be 
expected  as  the  oscillator  resonances  are  approached. 

For  completeness  we  show  In  Figures  9  and  10  the  results  of  methods  (3) 
and  (4),  respectively.  As  Is  seen,  the  variation  of  X^g^  with  composi¬ 
tion  X  Is  much  more  gradual  by  this  method  than  by  either  method  (2)  above  or 
the  approach  of  Section  4. 


d(An| 
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SECTION  6 


CONCLUSIONS 


Our  theoretical  oscillator  models  of  the  uniaxial  crystals  ZnS  and  CdS 
have  been  shown  to  yield  values  of  both  ordinary  and  extraordinary  refractive 
Indices  throughout  the  visible  and  near  Infrared  wavelength  regions  with  an 
accuracy  of  better  than  0.5Z.  However,  employing  the  theory  to  predict  the 
dispersive  properties  of  the  mixed  crystals  Zn^.^Cd^S  (l.e.,  the  Iso- 
Index  wavelength  and  slope  of  birefringence)  by  a  linear  Interpolation  has 
yielded  results  that  are  only  semlquantlatlvely  In  agreement  with  the  limited 
experimental  data  available. 

The  following  conclusions  characterize  the  results  of  our  study: 

*  The  Intricacies  of  band  structure  of  the  alloy  superlattices 
for  the  mixed  crystals  (Zn,Cd)S  are  too  subtle  to  allow  accur¬ 
ate  predictions  of  their  refractive  properties  by  simplified 
multiple  classical  oscllaltor  models. 

*  Better  results  are  obtained  when  e((o)  or  n^(«>),  rather  than 
n((»),  Is  represented  the  oscillator  model.  Predictions 
of  the  dependence  of  Iso-index  wavelength  and  slope  of  bire¬ 
fringence  upon  Zn  content  are  semlquantltatlvely  correct,  but 
only  so  long  as  we  Include  at  least  two  band-edge  oscillator 
resonances  In  our  mathematic  model. 

*  The  Fade  approximation  works  well  for  n((o),  provided  one  Inter¬ 
polates  directly  the  Fade  coefficients.  The  model,  however, 
yields  poles  coincident  with  the  Iso-lndex  point  at  unexpectedly 
low  levels  of  Zn  concentration. 

*  One  Indicator  that  shows  promise  Is  that  for  moderate  fractions 
of  Zn,  not  only  does  the  Iso-lndex  point  move  toward  shorter 
wavelengths,  but  according  to  our  approximate  theory  the  slope 
of  birefringence  Increases  above  that  of  pure  CdS.  This  sug¬ 
gests  the  possibility  of  more  highly  selective  Iso-lndex  fil¬ 
ters^  at  wavelengths  In  the  range  A600  A  to  4800  A. 

Finally,  we  believe  that  reliable  experiments  are  needed  to  yield  valid 
Information  on  the  desired  dispersive  characteristics  of  the  (Zn,Cd)S  alloys. 
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